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^ ■ ABSTRACT 

►~"S | We present an analysis of the long Rossi X-ray Timing Explorer (RXTE ) observation of the 

00 . Seyfert 1 galaxy MCG— 6-30-15, taken in July 1997. We have earlier used the data to place 

I ' constraints for the first time on the iron abundance - reflection fraction relationship, and now 

expand the analysis to investigate in detail the spectral X-ray variability of the object. Our 
CS| ■ results show that the behaviour is complicated. We find clear evidence from colour ratios 

and direct spectral fitting that changes to the intrinsic photon index are taking place. In gen- 
on : eral, spectral hardening is evident during periods of diminished intensity, and in particular, a 

general trend for harder spectra is seen in the period following the hardest RXTE flare. Flux- 
correlated studies further show that the 3—10 keV photon index T3_io steepens while that in 
the 10-20 keV band, rio-20, flattens with flux. The largest changes come from the spectral 
' index below 10 keV; however, changes in the intrinsic power law slope (shown by changes 

\ in r3_io), and reflection (shown by changes in rio-20) both contribute in varying degrees 

to the overall spectral variability. We find that the iron line flux F^a is consistent with being 
constant over large time intervals on the order of days (although the ASCA and RXTE spec- 
tra show that Fko changes on shorter time intervals of order < 10 ks), and equivalent width 
which anticorrelates with the continuum flux, and reflection fraction. A possible interpretation 
S-j \ for the iron line flux constancy and the relative Compton reflection increase with flux from 

c/3 . the flux-correlated data is an increasing ionization of the emitting disk surface, while spectral 

analysis of short time intervals surrounding flare events hint tentatively at observed spectral 
\ responses to the flare. We present a simple model for partial ionization where the bulk of the 

variability comes from within 6r g . Temporal analysis further provides evidence for possible 
' time (< 1000s) and phase (qb ~ 0.6 rad) lags. Finally, we report an apparent break in the 

power density spectrum (~ 4 — 5 x 10~ 6 Hz) and a possible 33 hour period. Estimates for the 
mass of the black hole in MCG— 6-30-15 are discussed in the context of spectral and temporal 
findings. 
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1 INTRODUCTION 

Observational and theoretical progress over the past decade have 
greatly improved our understanding of accreting black holes in 
Seyfert galaxies. Time-averaged X-ray spectra reveal a hard power- 
law continuum with a broad iron line and continuum reflection 
components (Pounds et al 1990; Nandra & Pounds 1994; Tanaka 
et al 1995; Nandra et al 1997). The power-law emission is pro- 
duced by thermal Comptonization in ~ 100 keV hot gas (Zdziarski 
et al 1994) above a thin accretion disc, which causes the reflection 

* Current Address: MIT Center for Space Research; 77 Massachusetts 
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(Guilbert & Rees 1988; Lightman & White 1988; George & Fabian 
1991; Matt, Perola & Piro 1991). 

Spectral variability in some objects show that the emission is 
produced in flarelike events which may move about over the in- 
ner accretion flow. The information from such variability is com- 
plex and its interpretation is not necessarily straightforward. Barr & 
Mushotzky (1986) and Wandel & Mushotzky (1986), by invoking 
the criterion for the fastest doubling time found in observations of 
AGN, demonstrate that a correlation between the X-ray luminosity 
and variability time scales exist and use this to place upper limits on 
the black hole masses in their sample of AGN observed with a slew 
of X-ray telescopes that include Ariel V, HEAO-1, and the Ein- 
stein observatory. Such a method however has been criticized for 
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its dependence on data quality and coverage (Lawrence et al. 1987; 
Mchardy & Czerny 1987). It does nevertheless provide a strong 
limit on the size of the emitting region. For example, Reynolds et 
al. (1995) and Yaqoob et al. (1997) have reported a factor of 1.5 
increase in flux in as little as 100s in MCG— 6-30-15. 

Others have resorted to employing power density spectrum 
(PDS) techniques (and variants thereof) as potential black hole 
mass estimators (e.g. Edelson & Nandra 1999 for NGC3516; 
Nowak & Chiang 1999 for MCG-6-30-15; Hayashida et al. 1998 
introduces the normalized power spectrum density (NPSD) for their 
sample of AGN observed with Ginga). However, even this is not 
without its caveats. X-ray variability studies of AGNs thus far (with 
the exception of IRAS18325-5926; Iwasawa et al. 1998) has shown 
that the power-law spectrum is essentially chaotic with no char- 
acteristic period. The problem of unevenly sampled data streams, 
characteristic of X-ray observations of AGN, is overcome by the 
solution of Lomb (1976) and Scargle (1982; see Press et al 1992). 

Others still have searched for time lags and leads using cross 
correlation function techniques. However, the problem associated 
with unevenly sampled data is still a concern, and has been ad- 
dressed by e.g. Edelson & Krolik (1988), Yaqoob et al (1997), 
Reynolds (1999), and in this paper. 

MCG-6-30-15 is a bright nearby (z=0.0078) Seyfert 1 galaxy 
that has been extensively studied by every major X-ray observa- 
tory since its identification. A recent study that takes advantage of 
the high energy and broad-band coverage of RXTE by Lee et al. 
(1998, 1999) with respectively a 50 ks observation in 1996, and 
400 ks 1997 observation have confirmed the clear existence of a 
broad iron line and reflection component in this object. A previous 
study by Iwasawa et al. (1996) (hereafter 196) using a 4.5 day ob- 
servation with ASCA revealed the iron line profile to be variable, 
which is confirmed in this long 1997 ASCA observation (Iwasawa 
et al. 1999, hereafter 199) which was observed simultaneously with 
RXTE . Additionally, Lee et al. (1999) have been able to constrain 
the iron abundance - reflection fraction relation for the first time 
using the RXTE observation. Guainazzi et al. (1999) give good 
bounds on the high energy cutoff of the continuum from a Bep- 
poSAX observation of MCG-6-30-15. 

We investigate in this paper changes in the direct and reflected 
components with our simultaneous RXTE and ASCA observations 
spanning a good time interval of ~ 400 ks, in order to speculate 
on possible causes for variability. In Section 3, we investigate rapid 
variability with colour ratios. This is followed in Section 4 by flux- 
correlated studies, and a detailed analysis of the time intervals sur- 
rounding the bright ASCA and RXTE flare shown in Fig. |IJ In Sec- 
tion 5, we shift gears to temporal studies and search (using cross- 
correlation techniques) for the presence of time lags and leads, in 
order to better understand the physical processes that may give rise 
to each other. We discuss in Section 6 tentative evidence for the 
presence of a possible 33 hour period, and break in the power spec- 
trum. Finally, we discuss the spectral and temporal findings in Sec- 
tion 7. In particular, we discuss the challenge that spectral findings 
may pose to current models for reflection by cold material, and im- 
plications from temporal studies for placing constraints on the size 
of the emitting region and mass of the black hole in MCG— 6-30- 
15. We present also a simple model to explain some of our enig- 
matic spectral findings. We conclude in Section 8 with a summary 
of the pertinent findings from this study. 



2 OBSERVATIONS 

MCG-6-30-15 was observed by RXTE over the period from 1997 
Aug 4 to 1997 Aug 12 by both the Proportional Counter Array 
(PCA) and High-Energy X-ray Timing Experiment (HEXTE) in- 
struments. We note that good data covered ~400 ks, even though 
the RXTE observation spanned ~700 ks. It was simultaneously ob- 
served for ~ 200 ks by the ASCA Solid-state Imaging Spectrom- 
eters (SIS) over the period 1997 August 3 to 1997 August 10 with 
a half-day gap part way through the observation. We concentrate 
only on the RXTE PCA observations in this paper. 



2.1 Data Analysis 

We extract PCA light curves and spectra from only the top Xenon 
layer using the FTOOLS 4.2 software. The light curves were ex- 
tracted with the default 16s time bins for Standard 2 PCA data 
time resolution. Data from PCUs 0, 1 , and 2 are combined to im- 
prove signal-to-noise at the expense of slightly blurring the spectral 
resolution. Data from the remaining PCUs (PCU 3 and 4) are ex- 
cluded because these instruments are known to periodically suffer 
discharge and are hence sometimes turned off. 

Good time intervals were selected to exclude any earth or 
South Atlantic Anomaly (SAA) passage occulations, and to ensure 
stable pointing. We also filter out electron contamination events. 

We generate background data using PCABACKEST V2.0C in 
order to estimate the internal background caused by interactions 
between the radiation/particles and the detector/spacecraft at the 
time of observation. This is done by matching the conditions of ob- 
servations with those in various model files. The model files used 
are the L7-240 background models which are intended to be spe- 
cialized for application to faint sources with count rate less than 
40ctss- 1 PCir 1 . 

The PCA response matrix for the RXTE data set was created 
using PCARSP V2.36. Background models and response matrices 
are representative of the most up-to-date PCA calibrations. 

Figure [j] shows the background subtracted light curve for 
ASCA and RXTE , with 256s binning. The time intervals corre- 
sponding to the ASCA and RXTE flares are decomposed further 
and analyzed in depth in Section 4.3. Significant variability can be 
seen in both light curves on short and long timescales. Flare and 
minima events are seen to correlate temporally in both light curves. 



3 RAPID VARIABILITY 

Intraday variability has been seen in most quasars and AGNs; in 
MCG— 6-30-15, rapid X-ray variability on the order of 100 s has 
been reported by a number of workers (e.g. Reynolds et al. 1995; 
Yaqoob et al. 1997). These significant luminosity changes on time 
scales as short as minutes can have strong implications for restrict- 
ing the size of the emitting region, and efficiency of the central en- 
gine. Colour ratios are good tools for discerning the processes that 
may give rise to these rapid variability effects. For the purposes 
of our study, we define the following count rate hardness ratios in 
order to assess the processes that may give rise to variability : 

5 - 7 keV 



Ri = 



R 2 = 



Ra = 



3- 


4.5 keV 


7.5 


- 10 keV 


3- 


- 4.5 keV 


10- 


- 20 keV 



3-A.5keV 



R 3 = 



Rs = 



7.5 - 10 keV 

5-7 keV 
10 - 20 keV 
5-7 keV 



(1) 



(2) 



(3) 



© 0000 RAS, MNRAS 000, 000-000 



The X-ray variability of the Seyfert 1 galaxy MCG—6-30- 1 5 from long ASCA and RXTE observations 3 

RXTE and ASCA light curves 



« 8 

■+-> 

d 
pi 
o 
o 

Eh ^ 

X 

K 



CD 



O 

o 
o 



C\2 




2x10 



4xicr 

Time (seconds) 



6x10" 



Figure 1. Background-subtracted light curve of MCG-6-30-15 for observations with 3 PCUs in the 2-60 keV band for the RXTE PCA, and 0.6-10 keV 
ASCA SIS. The epoch of the start and stop times for RXTE and ASCA respectively is 1997 August 4, 03:27:06 (UT) and 1997 August 12, 12:34:14 (UT) 
(RXTE ), and 1997 August 3 to 1997 August 10 (ASCA ). Light curves are binned using 256s time intervals. The intervals corresponding to the ASCA and 
RXTE flares are further subdivided and discussed in depth in Section 4.3. 



Hardness ratio 



Rl 


3.7 


R2 


3.25 


R3 


0.95 


R4 


4.42 


R5 


1.95 


R6 


1.24 



Table 1. Assessment for amount/absence of variability for the hardness ra- 
tios shown on the left panels of Figs. ^|— ^. x 2 , is evaluated for 138 degrees 
of freedom. 



Re 



10 - 20 keV 
7.5 - 10 keV 



We expect the reflection component to dominate above 10 keV. At 
the lower energies, the dominant features include the iron line be- 
tween 5.5-6.5 keV and the warm absorber below 2 keV. We note 
also that the iron line contribution to the overall flux in the 5-7 keV 
band is only ~ 15%, and hence be hereafter referred to as the iron 
line band. 

The Ri assessment of the (5-7 keV) iron line band versus (3- 
4.5 keV) lower continuum (Fig. ^a) shows that the source is in- 



trinsically harder during the minima. Similarly, R2 hardness ra- 
tio of the (7.5-10 keV) upper continuum to the (3-4.5 keV) lower 
continuum reveals a similar trend (Fig. ^a); likewise, R4 assess- 
ment of the (10-20 keV) reflection hump and (3-4.5 keV) lower 
continuum (Fig. mi). R3, R5, R6 show no obvious trends (respec- 
tively Figs. Bp,kpMc) from count-rate versus hardness ratio plots, 
although subtle correlations can be seen in the time versus hardness 
ratio plots. In particular Figs, ^b and ^c show dramatic hardening 
of the spectrum during the time periods following the RXTE flare. 
(We mark the beginning of the ASCA and RXTE flares respectively 
with 'A' and 'X' in the light curves shown in Figs. ^.) In order 
to further quantify the amount / absence of variablity, we apply a 
X 2 test to the hardness ratio trends as shown in the left panels of 
Figs.^— Table [l] details these results for 138 degrees of freedom. 



3.1 The unusual properties of the deep minimum period 
following the hard RXTE flare 

We mark the regions corresponding to the ASCA and RXTE flare 
events in Figs. ^ ^. The time span corresponding to the ASCA and 
RXTE flares were chosen in order that the times surrounding the 
flares and minimum (following said flare) can be compared. The 
spectral behaviour during these periods is thoroughly investigated 
in Section 4.3. For present purposes, we draw the reader's atten- 
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tion to the deep minimum immediately following the bright RXTE 
flare. (The general location of this minimum in the light curve is 
labeled with 'DM' in Figs. ^.) Despite the hardness ratio re- 
sults discussed previously (e.g. no obvious trends seen in R3, R5, 
Re), the spectrum during this time period always hardens in all 
colours, whereas this is not necessarily the case with the other min- 
ima events (e.g. the deep minimum following the soft ASCA flare). 
In addition, it is not clear whether the spectrum reaches its lowest 
minimum slightly before or after an observed hardening. 

We also bring attention (by means of a horizontal line through 
the hardness ratio versus time plots in Figs. ^— §) to the general 
trend for harder spectra in the ~ 260 ks time period that begins 
with the RXTE flare. 



3.2 Possible causes of Hardness Ratio Variability 

Many complicated processes contribute to the overall temporal and 
spectral behaviour of an AGN. It is necessary to disentangle these 
components in order to assess the physical processes that are re- 
sponsible for the observed variability phenomenon. Effects due to 
reprocessing are significant only at the higher energies, the (10- 
20 keV) band being the most sensitive probe of any subtle changes 
in the reflection component. Martocchia & Matt (1996) suggested 
that gravitational bending/focusing as the X-ray source gets closer 
to the black hole will enhance the amount of reflection. Therefore, 
it is possible that flares closer to the hole during the minima can en- 
hance the amount of reflection, through increased beaming of the 
emission towards the disk. We conclude however from the com- 
bined hardness ratio findings described thus far that the spectral 
hardening during periods of diminished intensity point largely to 
changes in the intrinsic photon index T being the main culprit for 
the observed spectral variability. We are led to this conclusion be- 
cause spectral changes are seen in all bands, which is most likely 
to be due to changes in the spectral index. In other words, we may 
be seeing the effects of changes in the coronal temperature affect- 
ing the intrinsic power law slope, or that changes in F are due to 
coverage of flares such that we are observing the effects of flares 
occurring at different heights (e.g. Poutanen & Fabian 1999). How- 
ever, as we shall now show, direct spectral analysis reveals a more 
complex situation. 



4 SPECTRAL FITTING 

Having established that variability effects are many, we next inves- 
tigate spectral changes in time sequence in order to obtain a better 
understanding of variability phenomena between the different flux 
states with particular emphasis on the flares and deep minima. Data 
analysis is restricted to the 3 to 20 ke V PCA energy band. The lower 
energy restriction at 3 keV is selected in order that the necessity for 
modelling photoelectric absorption due to Galactic ISM material, 
or the warm absorber that are known to be present in this object is 
removed. (Lee et al. 1999 have shown that effects due to the warm 
absorber at 3 keV are negligible for this data set.) 

4.1 Spectra selected in time sequence and by flux 

We choose the time intervals (Fig. |l|) in order to contrast and study 
the different variability states. They are defined such that intervals 

11 and i6 correspond to the two RXTE deep minima, and intervals 

12 and i7 for the relatively calm periods following these minima, 



and i5 a flare event for comparison. We investigate in detail in Sec- 
tion 4.3 the times surrounding the RXTE counterpart of the soft 
ASCA flare, and the hard flare observed by RXTE (unfortunately 
missed by ASCA ), as depicted in Fig. [j] and subsequently in Figs|^, 
and [u|. Both the ASCA and RXTE light curves for the full obser- 
vation are shown in Fig. [j]. 

We also separate the data according to flux in order to assess 
whether a clear picture can be developed of the dominant processes 
that may be at work for a given flux state. To do so, we separate the 
400 ks observation by flux, with f2 and f3 being the intermediate 
states between the lowest fl (2.84 x 10 -11 ergcm~ 2 s _1 ), and 
highest f4 (4.79 x 10 -11 erg cm -2 s _1 ) fluxes. Properties of these 
flux levels (in the 2-60 keV, and 3-20 keV energy bands), and il i7 
intervals are detailed respectively in Table 0, and Table IvL 



4.2 Temporal Variability and Spectral Features 

A nominal fit to the entire data set (ie. ASCA , RXTE PCA and 
HEXTE) demonstrated the clear existence of a redshifted broad 
iron Ka line at ~ 6.0 keV and reflection hump between 20 and 
30 keV as shown in Lee et al. (1998, 1999). In general, fits to the 
solar abundance models of George & Fabian 1991 (hereafter GF91) 
are better than simple power law fits and further reinforce the pref- 
erence for a reflection component. Since the reflected continuum 
does not contribute significant flux to the observed spectrum be- 
low 10 keV, a simple power law and iron line fit to the data below 
10 keV will reveal changes in the intrinsic power-law of the X- 
ray source. If indeed spectral variability is due to changes in the 
intrinsic power law slope which would implicate changes in the 
conditions of the X-ray emitting corona, we should see noticeable 
changes in the values of r( 3 _ 10 ) between the different temporal 
states. However, because reflection only significantly affects ener- 
gies above 10 keV, changes seen in r( 10 _2o) would suggest that 
variability may be due to the amount of reprocessing. (We have 
tested this premise by comparing simulated spectra in XSPEC us- 
ing the PEXRAV model and find that the overall flux contribution 
from reflection alone is > 60 per cent above 10 keV.) This would 
have strong implications for geometry (ie. where the direct X-ray 
flares are partially obscured), motion of the source (eg. Reynolds & 
Fabian 1997; Beloborodov 1998), or gravity (light-bending effects 
that will beam/focus more of the emission down towards the disk; 
Martocchia & Matt 1996). 



4.2.1 The different flux states 

We first investigate whether a correlation exists between flux and 
the various fit parameters, and in particular whether changes in re- 
flection are dramatic. We do so by first fitting a simple model that 
consists of a power law and redshifted Gaussian to the 3-10 keV 
band, and comparing that with a power law fit to the 10-20 keV 
band. In doing so, we find that while the intrinsic 3-10 keV power 
law slope increases, Tio-20 appears to flatten with increasing flux 
(Fig. ^). Additionally, the iron Ka flux Fk« does not change with 
any statistical significance while the flux between the lowest fl and 
highest f4 flux nearly doubles (Table ^|); similar results for the con- 
stancy of the iron line in MCG— 6-30-15 was noted by McHardy 
et al. (1998), and by Chiang et al. (1999) for Seyfert 1 galaxy 
NGC5548. The ratio plot of data-to-model using a simple power 
law fit to the 3-20 keV data clearly illustrates the difference in the 
line and reflection component between the lowest and highest flux 
states in Fig. Ra. 
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Figure 3. Hardness ratio plots against time and count rate of the upper continuum (7.5-10 keV) versus (a) lower continuum (3-4.5 keV) and (b) iron line band 
(5-7 keV). Light curves are binned into 5000s intervals, and represent the sum of the light curve of the two energy bands being compared.. 
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Figure 4. Hardness ratio plotted against time and count rate for the reflection hump (10-20 keV) versus (a) the lower continuum (3-4.5 keV), (b) iron line 
band (5-7 keV), and (c) upper continuum (7.5-10 keV). Light curves are binned into 5000s intervals, and represent the sum of the light curve of the two energy 
bands being compared.. 
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Change in r for different states of MCG-6-30-15 



Flux 
Levels 


pa 

1 3-10 


F b 

Ko 

ph cm -2 s -1 


1 10-20 


Ka 

eV 


Flux e 
(3-10 keV) 


Exposure 
10 4 s 


2-60 keV 

ct s _1 


Flux cuts 

ct s -1 


Time range 
10 4 s 


fl 
£2 

a 

f4 


1.91 ±0.02 
1.93 ±0.02 
1.95 ±0.01 
1.98 ±0.01 


-1 H1 +0.16 
1 89 +0 - 19 

i 79 +0.14 
1 -' z -0.14 
i ocr + 0.19 

l.Bt)_ 17 


1.99 ±0.10 
1.65 ±0.09 
1.55 ±0.07 
1.48 ±0.07 


466+ 41 

wu -36 
qq.0,+39 

272J_ 2 J 


2.84 
3.44 
3.94 
4.79 


7.8 
5.7 
7.7 
6.2 


9.89 ± 0.34 
13.03 ± 0.41 
15.43 ±0.36 
19.33 ± 0.40 


f < 12 
12 < f < 14 
14 < f < 17 

f > 17 


full 
full 
full 
full 



Table 2. Results are quoted from simple power law fits. a Power-law photon index for 3 keV < E < 10 keV. b Flux of iron emission line in units of 10 — 4 
phcm - 2 s -1 . c Power-law photon index for 10 keV < E < 20 keV. d Equivalent width of the iron emission line. e 3-10 keV flux in units of 10 — 11 
erg cm -2 s — 1 Intervals fl-f4 correspond to the full time-averaged spectrum separated according to flux. Count rates are given for the 2-60 keV energy band 
offheftXre PCA. 



Spectral Fits using power-law with reflection model for different flux states of MCG— 6-30-15 
Data a r 4 _20 b A c refl d LineE c a * F Ka §W (eV) h Flux l \ 2 



fl 

f2 

a 

f4 



1.88 
1.97 
2.01 
2.07 



0.03 
0.03 
0.03 
0.03 
+0.03 
0.02 
+0.03 
0.02 



1.16 
1.57 
1.91 
2.51 



0.03 
0.04 
+0.05 
0.05 
+0.05 
0.05 
+0.09 
0.05 



0.35 
0.91 
1.10 
1.37 



+0.15 

0.16 
+0.22 
0.20 
+ 0.19 
0.16 
+0.23 
0.14 



6.04 
6.03 
6.05 
5.98 



+0.05 
0.05 
+0.07 
0.07 
+0.06 
0.06 
+0.09 
0.10 



0.54^ 



54+ ' 07 

u - oy: -0.06 
4+0.09 
-0.10 

u - oa -o.io 
36+° 14 



1.50 
1.31 
0.96 
0.72 



+0.11 
-0.12 

0.16 

0.16 
+0.12 

0.11 
+0.14 

0.19 



375: 
26i; 
167: 
98+ 



-30 
-31 
-31 
-21 
-1!) 
19 



4.61 
5.73 
6.58 
8.03 



76 
44 
34 
35 



Table 3. a Power-law photon index. b Power-law flux at 1 keV, in units of 10 — 3 phcm - 2 s -1 keV - , c Reflective fraction = fl/2ir. d Energy of the iron 
Ka emission line. e Line width in units of eV. * Flux of iron emission line in units of 10~ 4 phcm -2 s _1 . 9 Equivalent width of the emission line. h 
3-20 keV flux in units of 10 — 11 erg cm -2 s~ 1 . 1 \ 2 f° r 39 degrees of freedom. 



Having demonstrated the existence of a strong reflection com- 
ponent in Fig. ^a, and assessed its nature with more complicated 
fits in Lee et al. (1999) for this data set, we next investigate the fea- 
tures of reflection in detail for the different fluxes by fitting the data 
with a multicomponent model that includes the reflected spectrum. 
The underlying continuum is fit with the model PEXRAV which is a 
power law with an exponential cut off at high energies reflected by 
an optically thick slab of neutral material (Magdziarz & Zdziarski 
1995). We fix the inclination angle of the reflector at 30° so as to 
agree with the disk inclination one obtains when fitting accretion 
disk models to the iron line profile as seen by ASCA (Tanaka et al. 
1995). Due to the strong coupling between the fit parameters of T, 
abundances, and reflection, we fix the low-Z and iron abundance 
respectively at 0.5 and 2 solar abundances as determined by Lee et 
al. (1999) for the fits presented in Table [| the high energy cutoff is 
fixed at 100 keV appropriate for this object (Guainazzi et al. 1999; 
Lee et al. 1999). An additional Gaussian component is added to 
model the iron line. 

Using this complex model, we find that the 4-20 keV power 
law slope and reflection fraction R increases with flux (Table [|) 
while the strength of the iron line, Jka decreases (Fig. ^i), the lat- 
ter in contrast to the findings for a constant fka discussed in the 
context of simpler fits. (Fkc is defined as the total number of pho- 
ton flux in the line.) We note that Fko is consistent with constancy 
if unity abundances are assumed; this is in agreement with simple 
power law fits. However, this leaves us with difficult-to-constrain 
errors, and worse fits in a \ 2 sense. (With the exception of Fkc, all 
other parameters as e.g. F shown in Table |3| follow similar trends 
whether unity or non-unity abundances are assumed.) It is clear that 
degeneracies exist and cannot be resolved with 3—20 keV RXTE 
data - we suspect that the dependence on abundance is largely due 
to the modelling of the iron edge in these complex fits. Neverthe- 




Figure 5. Intrinsic photon index T3_io vs Tiq— 20 of reflection component 
shows that changes in both are evident. 



less, we test this hypothesis by including an edge feature to the sim- 
ple power law model of Table ^. These results presented in Table ^ 
show that FKa is indeed consistent with constancy and strength- 
ens the argument that the behaviour of the iron line as given by the 
complex model of Table |3| is complicated with degeneracies. This 
and the flux constancy of the iron line is well illustrated in Fig. 
which shows the ratio of the best-fit data against the model of Ta- 
ble j| Certainly, the case is strong for a requirement of supersolar 
iron abundances and is reflected in the strength of the iron line. We 
discuss this in depth in Lee et al. (1999). 

While FKa apparently decreases with flux from complex fits, 
we find that the reflection fraction, and absolute normalization of 
the reflection component (Rnorm) increases with flux (Table ^ 
and Fig^b). We define Rnorm — A * R, where A is the power 



© 0000 RAS, MNRAS 000, 000-000 



8 J. C. Lee et al. 




(b) 



(c) 




2x10 ° 3x10 
Absolute Reflection 



(d) 



r (3-10) 



0.5 1 

Reflection Fraction 



Figure 6. Plots representing results from complex fits using the model : PEXRAV + GAUSSIAN discussed in Section 4.2.1, and detailed in Table ^[ 'Absolute 
reflection' refers to the absolute normalization of the reflection component (in units of 10 -2 phcm -2 s _1 keV -1 ). 
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Figure 7. (a) Ratio plot of data-to-powerlaw_model of the iron line and reflection component for comparison between the extreme flux states fl (lowest) and 
f4 (highest). A simple power law is used to fit the data, (b) Ratio of bestfit data-to-model for the flux states fl and f4 illustrates the flux constancy of the iron 
line. The model used is that presented in Table |^. Dotted lines represent fl and solid lines f4. 
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Spectral Fits using power-law with edge model for different flux states of MCG— 6-30- 1 5 
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Table 4. a Power-law photon index. 6 Energy of the iron Ka emission line. c Line width in units of keV. d Flux of iron emission line in units of 10~ 4 
phcm~ 2 s _1 . e Equivalent width of the emission line (eV). f Energy of absorption edge (redshifted to z=0.()078 appropriate for MCG— 6-30-15) in units of 
keV. 9 Maximum value for the optical depth at threhold energy. h x 2 for 12 degrees of freedom. 



law flux at 1 keV, in units of 10~ 3 phcm" 2 s" 1 keV -1 . (Fig. ^a 
clearly shows that stronger reflection is present during higher flux 
states.) The anticorrelation between _Fka and reflection (Fig. ^c) 
can be due to an 'artificial' effect in which the presence of a 
strong reflection spectrum during the high flux states has the ef- 
fect of removing part of the flux from the line in the fits, thereby 
resulting in lower observed Fna during the higher flux states. 
This is linked to the strong coupling between the fit parameters of 
F, R, and abundance, discussed previously. There is the possibil- 
ity that iron becomes more ionized as the flux increases which will 
weaken the observed line flux. (We discuss ionization scenarios in 
Section 7.) As always we caution the effect of inadequate spectral 
resolution and possibly incomplete models. 

Additionally, we find that Wko anticorrelates with R 
(Fig. ^d). This lack of proportionality between R and Wkc is un- 
expected in the context of the standard corona/disk geometry, the 
implications of which are discussed in Section 7. Chiang et al. 
(1999) find similar results in their multi- wavelength campaign of 
NGC5548. In light of degeneracies associated with complex fits, 
we primarily present our results from simple power law fits. 

4.3 The bright ASCA and RXTE flares 

We next investigate in greater detail the time sequences surround- 
ing the brightest ASCA and RXTE flares. In particular, we are 
prompted by the peculiar behaviour of the iron line surrounding 
the time interval of the ASCA flare as reported by 199 for this time 
sequence as seen in the ASCA data. According to 199, there is dra- 
matic change of the line in both profile and intensity during this 
period: interval a as depicted in Fig. |s| is marked by an extremely 
redshifted line profile that is characterized by a sharp decline in the 
line energy at 5.6 keV (far below the 6.4 keV rest energy of the line 
emission) with a red wing that extends to 3.5 keV, and line inten- 
sity that is ~3 times that for the time averaged data. The transition 
to time interval b is marked by an abrupt factor of 2.2 drop in the 
averaged 0.6-10 keV count rate, and similar ~2 drop in intensity. 
(We present in this section the behaviour as seen by RXTE .) Fig. H 
shows the ASCA and corresponding RXTE light curves for this 
period of interest. We note that the two observations are slightly 
offset from each other; this will allow for a good assessment of 
the events immediately preceding and following this flare event. 
Unfortunately ASCA data do not exist to coincide with the most 
prominent RXTE flare during the time interval between 2.8 x 10 
and 3.8 x 10 5 s. 

We fit the 3-10 keV data with a number of different mod- 
els : (1) power law + Gaussian, (2) (power law + Gaus- 
sian) modified by absorption, (3) POWER LAW + DISKLINE and 
(4) (POWER LAW + DISKLINE) modified by absorption. For the 



majority of the fits, the DISKLINE model by Fabian et al. (1989) 
provided the best results; these POWER LAW + DISKLINE best fit 
values are detailed in Table | and Table [], respectively for the ASCA 
and RXTE flare events. (We note however that differences between 
X 2 for fits to the different models are not great : Ax 2 less than 2 
for 1 extra parameter.) Due to the inadequate resolution of RXTE 
in the iron line band such that it is insensitive to the details that a 
diskline model can provide, we fix DISKLINE parameters at best fit 
ASCA values reported by 199 for this 1997 observation : emissivity 
a — —4.1, inclination i — 30°, respectively inner and outer radius 
Rin = 6.7, and R ou t = 24. We caution that fixing these values 
is likely to be an oversimplification of the true scenario since the 
line profile is known to vary on short time scales (e.g. Iwasawa et 
al. 1996), the impact of which is not easily assessable given present 
data quality. Nevertheless, we find that it is presently the best op- 
tion for giving a first order approximation of what is likely to be 
happening. 



4.3.1 The ASCA flare 

Figs. |i] & ^ illustrate the RXTE time intervals corresponding to the 
periods surrounding the ASCA flare, with periods Al and A2 cho- 
sen to respectively correspond to the times immediately preceding 
and following the flare 'a' seen in ASCA . Table ^ details fit results 
for the time intervals depicted in Fig. |^. 

There are noticeable changes in the intrinsic power law slope 
throughout the intervals of interest. In particular, the most dra- 
matic changes occur between the intervals q and e. T3_io steep- 
ens from q to Al intervals immediately preceding the flare event 
seen in ASCA , which may have the effect of producing the no- 
ticeable decrease in Wkc ■ In the ~ 1 1 ks separating the ' pre- ' 
and ' post- ' flare event (respectively Al and A2), Fkc increases 
(from 1.67 ± 0.44 to 2.61 ± 0.38 phcm" 2 s" 1 ) and Wk<* (from 
278 ± 74 to 543 ± 79 eV), nearly double. This is consistent with 
199 findings from ASCA data for factor of ~ 2 difference in Fkc 
for flare a, and post-flare b intervals shown in Fig. |§| It also appears 
that Fio-20 steepens during this transition, although errors are too 
large to definitively make this claim. Subsequent changes to T3_io 
are still noticeable ; Fio-20 appears to flatten in e, and along with 
Wku and flux remain fairly constant through g. Unfortunately, the 
errors associated with Fko are generally too large to make statisti- 
cally significant statements about changes in the iron line flux, even 
though there are indications for this (e.g. Fig. [kJ; notice especially 
the difference in Fko between Al and A2. 
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Change in r for time periods previous to and following ASCA flare 



Interval 


pa 

1 3-10 


Ka 


F c 

Ka 


pd 

1 10-20 


W Ka (eV) 


Flux f 


q 

Al 


1.92 ± 0.04 
2.05 ± 0.04 


6.50 ±0.15 
6.37 ±0.19 


1.88 ± 0.35 
1.67 ±0.44 


, c-a+0.25 
1 - OZ -0.27 
, „+0.24 
1 - do -0 .27 


466 ± 87 
278 ± 74 


3.25 
4.62 


A2 


1.99 ± 0.04 


6.38 ±0.11 


2.61 ± 0.38 


1 ni+0.28 
i - ai -0.27 


543 ± 79 


3.80 


e 


1.85 ± 0.04 


6.45 ±0.11 


2.26 ± 0.32 


1 cc+0.24 
L66 -0.27 


566 ± 80 


3.18 


f 


1.97 ±0.03 


6.47 ±0.09 


2.34 ±0.27 


1 70+ - 18 


471 ± 55 


3.98 


g 


1.96 ± 0.03 


6.48 ±0.10 


2.42 ± 0.30 


1 61+ 019 
^•Di_0.18 


466 ± 58 


4.17 



Table 5. Results are quoted from simple power law fits; the diskline model of Fabian et al. (1989) is used to model the line emission. a Power-law photon 
index for 3 keV < E < 10 keV. b Centroid energy of the iron Ka line. c Flux of iron emission line in units of 10 4 phcm 2 s 1 . d Power-law photon index 
for 10 keV < E < 20 keV. d Equivalent width of the iron emission line. e 3-10 keV flux in units of 10 — 11 erg cm -2 s" 1 . 
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Figure 8. RXTE (top) and ASCA (bottom) light curves corresponding to the 
most prominent ASCA flare. Intervals a and b in bottom panel correspond 
to the flare and post-flare phases discussed in Iwasawa et al. (1999). 
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Figure 9. RXTE light curve depicting the intervals of interest preceding 
and following the most prominent flare a seen by ASCA . 



4.3.2 The RXTE flare 

We next investigate the time intervals surrounding the RXTE flare 
(Fig. |ll]); due unfortunately to an absence of ASCA data during 
this time interval, we cannot make analogous comparisons. How- 
ever, we find that changes to the intrinsic power law slope during 
these intervals follow similar trends presented for the RXTE coun- 
terpart of the ASCA flare in the previous section, although events 
surrounding this bright RXTE flare event appear to be much more 
erratic and complicated. (We remind the reader that this is also ap- 
parent in hardness ratio comparisons during this time interval.) The 
similarities with the ASCA flare lie in observed trends to changes 
in F3_io. In particular, there is a noticeable flattening of F3_io in 
the transition between ' pre- ' and ' post- ' flare states XI to X2, 
which continues through interval m, before suddenly steepening in 
the following interval n. Unfortunately, errors are such that again 
we are unable to make any statements regarding changes to Jka 
or Tio-20- However, we note that Fkc* is generally high for these 
time intervals associated with this hard RXTE flare event. Table ^ 



F (10 11 erg cm 3 s 1 ) 

Figure 10. RXTE light curve depicting the intervals of interest preceding 
and following the most prominent flare seen by ASCA . 
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Figure 11. RXTE light curve depicting the intervals of interest preceding 
and following the most prominent flare seen by RXTE . 
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Figure 12. RXTE light curve depicting the intervals of interest preceding 
and following the most prominent flare seen by RXTE . 



details these results. As with the ASCA flare, Fig. [l§for the RXTE 
flare hints at changes in fk Q on short time intervals (e.g. XI and 
m) 



4.4 The intervals il to i7 

It is clear from the spectral analysis thus far that conditions can alter 
suddenly and erratically. In order to assess whether a more simpli- 
fied picture exists, we investigate spectral features of the deep min- 
ima in contrast with flare type events, using a model that consist 
of simple power law and redshifted Gaussian component. Table || 
confirms the findings of Section 4.2.1 such that in general, T^-io 
tends to be flatter during the minima in contrast to the flare states. 
A close comparison of T3_io versus Tio-20 for the differing states 
suggests that we are largely seeing intrinsic changes in the power 



Interval 


Exposure 


2-60 keV 


Time range 




10 3 s 


ct s" 1 


10 4 s 


il 


6.74 


9.59 ± 0.12 


3.0-4.8 


i2 


9.52 


15.86 ±0.10 


6.0-8.0 


i5 


5.50 


15.78 ±0.14 


51.0-53.5 


i6 


17.65 


10.50 ± 0.07 


54.0-58.0 


i7 


19.95 


16.47 ± 0.07 


62.0-66.0 



Table 7. il-i7 correspond to the nine intervals of interest for the light curve 
given in Fig. |lj The intervals il and i6 correspond to the deep minima, 
whereas i2, i5, and 17 have comparable fluxes. Count rates for the intervals 
are given for the 2-60 keV energy band of the RXTE PCA. 



law slope rather than reflection, although it is likely that we are see- 
ing contributions from both effects. Additionally, ratio plots of data 
against model using a power law fit show that there is a noticeable 
change in the line flux, profile, as well as the reflection component, 
similar to that seen in Fig. ^a. 



4.5 Summary of spectral findings 

We find evidence from flux-correlated studies that r3_io steepens 
significantly with flux (Ar3_io ~ 0.06 for a doubling of the flux 
from fl to f4; Fig. |l3|a) while surprisingly, the iron line strength 
appears to remain constant (at most differing by ~ 1.5 x 10~ 5 
phcm -2 s _1 from flux-correlated studies). Changes to Tio-20 
(Ario-20 ~ 0.3) and W Ka (AW Ka ~ 200 eV) are also evident 
and anticorrelate with flux (Fig. |l3| b for the latter). 

A close look at events corresponding to deep minima versus 
flares reinforces the finding that changes to the intrinsic power law 
slope are evident with a comparably steeper Ta-io value during 
the flares. Figs. 1131 illustrate that the behaviour of the intrinsic pho- 
ton index and Wkc during the flares is consistent with the flux- 
correlated behaviour. 

We find that reflection increases with flux when fitting 
flux-separated data with a complex model that includes the re- 
flected spectrum. Curiously, reflection fraction R anticorrelates 
with Wife*; similarly the absolutely normalization of the reflec- 
tion component (Rnorm) anticorrelates with Jka- We note that 
contrary to findings from simple power law fits, fka is seen to 
decrease with flux, when non-unity abundances are assumed. We 
caution however about the large degeneracies associated with com- 
plex fits when only the 3—20 keV RXTE data are considered. (The 
case for requiring supersolar abundances is discussed in Lee et al. 
1999). With the exception of the behaviour of Fko, all other pa- 
rameters as e.g. F display similar trends using complex and simple 
power law fits. The apparent discrepancy between fka from sim- 
ple and complex fits in conjunction with the large errors associated 
with Fkci would suggest that an ambiguity exists in measurements 
of the iron line itself, which RXTE is unable to resolve. 

A detailed investigation of the time intervals surrounding the 
bright ASCA and RXTE flare reveal further complexities. While 
changes to r3_io are consistent with flux-correlated studies, we 
find that there is evidence to suggest that a change in Fko oc- 
curs during the time intervals immediately associated with the flare 
event; in other words, tentative evidence for changes to Fna are 
apparent on short time scales. In particular, a large increase in the 
line flux (e.g A2) is evident in the interval immediately following 
the flare; Wna increases similarly. It is curious that fka shows a 
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Change in r for time periods previous to and following RXTE flare 
Interval E b FSi rf n 9n iy« (eV) Flux f 

3—10 Ka Kq 10 — 20 Kq v ' 

v 1.97 ±0.03 6.47 ±0.14 2.13 ±0.37 1.74±°H 361 ± 63 4.65 

XI 2.05 ±0.03 6.40 ±0.12 2.50 ±0.38 1-73^;^ 357 ± 54 5.46 

X2 1.93 ±0.03 6.35 ±0.10 2.19 ±0.30 1.69±g;|| 409 ± 57 4.09 

m 1.83 ±0.03 6.49 ±0.10 1.94 ±0.24 l-57±°;^ 457 ± 56 3.35 

n 1.98 ±0.03 6.46 ±0.11 2.22 ± 0.30 1-83^°;^ 497 ± 68 3.58 

p 1.94 ±0.02 6.45 ±0.07 2.48 ± 0.22 l.Blto^* 501 ± 45 3.96 

Table 6. Results are quoted from simple power law fits; the diskline model of Fabian et al. (1989) is used to model the line emission. a Power-law photon 
index for 3 keV < E < 10 keV. b Centroid energy of iron Ka line. c Flux of iron emission line in units of 10 4 phcm 2 s 1 . d Power-law photon index 
for 10 keV < E < 20 keV. e Equivalent width of the iron emission line, 3-10 keV flux in units of 10 — 11 ergcm~ 2 s _1 . 



Change in r for different states of MCG-6-30-15 
FL r io-20 W d ( e V) Flux" 

11 1.80 ±0.06 l.79to.4B 1-77 ±0.30 5391^ 2.47 

12 1.97 ± 0.03 1.73±q;29 1.59 ±0.19 322+H 3.98 

15 1.94 ±0.03 1.91+2;^ 1.58 ±0.16 318t^ 3 4.17 

16 1.82 ±0.03 1.85tn|| 1.56 ±0.16 494^ 2.68 

17 1.98 ±0.02 1.98+Q 32 1-83 ±0.14 305+gg 4.33 

Table 8. Results are quoted from simple power law fits; a redshifted Gaussian is used to model the line emission. a Power-law photon index for 
3 keV < E < 10 keV. b Flux of iron emission line in units of 10~ 4 phcm" 2 s" 1 . c Power-law photon index for 10 keV < E < 20 keV. d Equivalent 



significant increase after the flare rather than during, and may be an 
indication that are we witnessing some type of response to the flare 
(e.g. Fig. |l(^), but also caution that the evidence is very tentative. 
We note also that Fk is comparably high (~ factor of 1 .7 increase) 
during the times surrounding the RXTE flare, and times following 
the ASCA flare, in contrast to values presented in Tables ^, and 

1 



5 TIME LAGS, LEADS, AND REVERBERATION 

Motivated by the temporal findings of e.g. Miyamoto et al. (1988) 
and Cui et al. (1997) for Cygnus X-l, we next wish to investigate 
whether the collecting area of RXTE coupled with this long obser- 
vation is sufficient to discern time lags, time leads, and in particular 
whether reverberation effects can finally be seen. We acknowledge 
that a good assessment can be hampered by unevenly sampled data, 
the nature of which has been investigated by a number of work- 
ers (e.g. Edelson & Krolik 1998; Yaqoob et al. 1997) for AGN to 
In't Zand & Fenimore (1996) for application to gamma-ray bursts. 
We adopt a method similar to that presented by Edelson & Kro- 
lik (1998). Accordingly, we define the following formalism for our 
calculations of the autocorrelation function (ACF), and cross cor- 
relation function (CCF). 

CCF(n) = ±-Y y ' ii)Z ' {i + T) for (4) 

i— — n 

V"_ z . 

where z = — — 

n 



and CCF(T ) = ^2(yi-y){zi-z) = M 

i=l 

The time interval r is incremented by n multiples of 64s (r = 
n x 64); the subsequent definition of n T is the number of bins 
for which the difference in time between consecutive time intervals 
satisfy the present value for r. The variables K and M correspond 
respectively to the ACF and CCF value at r = 0; these terms are 
used in order to normalize the ACF and CCF so that coherent noise 
addition at r = is eliminated. 

In order to better understand the nature of our findings shown 
in Figs. [T^ - |l7| , we run our CCF algorithm on simulated light curves. 
The light curves over which the CCFs are evaluated are identical 
except that one has a specified fraction shifted in phase or time. In 
principle, time- and phase-shifted light curves are identical at cer- 
tain Fourier frequencies as for example, for the case in which the 
power spectrum can be represented by a single sine curve. How- 
ever, if the light curve is the addition of e.g. several sine curves 
(more representative of actual physical systems), then a time shift 
will have the effect of shifting the entire pattern by the specified 
time interval, whereas a phase shift will shift each individual sine 
curve along by that phase. The final outcome of the phase-shifted 
light curve will consist of the additive components of the individ- 
ual sine curves that have been shifted, i.e. the individual frequen- 
cies are added together separately. The CCF of the simulated phase- 
and time-shifted light curves are shown in Figs.[H|; solid lines cor- 
respond to the CCF, and dashed lines to the mirror image of this. 
We point out that the CCFs show subtle differences. For the phase- 
shifted light curve of Fig. [l4|a, a constant offset from zero is seen 
between the CCF and its mirror image (represented by respectively 
the solid and dashed lines). This is contrasted with the time-shifted 
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Figure 13. 2—10 keV flux versus (a) T3_io, and (b) equivalent width Wkq- Filled circles represent the different flux states; filled squares and dashed lines 
represent the ASCA flare intervals; triangles and dot-dashed lines represent RXTE flare intervals. 



light curve in which the CCF and its mirror image become compa- 
rable (i.e. there is no persistent offset) at ~ 30 bins (i.e. ~ 1900s). 
Simulated light curves are generated using Monte Carlo techniques 
for the flux, with power inversely proportional to frequency, while 
times are identical to those in the real data. We note that the simi- 
larity between the phase- and time-shifted data between and ~10 
bins (i.e. ~ 0-640 s) is due to the fact that only 20 per cent of the 
light curve has been shifted, with the other 80 per cent of the light 
curve that remain identical between the phase-shifted and time- 
shifted simulations. 

Assessment of the CCFs between the upper continuum (E3 : 
7.5-10 keV) and lower continuum (El : 3-4.5 keV) and iron line 
region (E2 : 5-7 keV) shown in Fig. [^, reveal evidence for a pos- 
sible phase shift. In a comparison with the CCF of the simulated 
light curve shown in Fig. |l4p , in which an artificial phase lag of 
<f> ~ 0.6 rad is introduced, we find that a similar trend exists in the 
actual data. This would suggest that ~ 10-20 per cent of the upper 
continuum band lags that of the lower energy bands El and E2. 

We next assess the nature of the (10-20 keV) reflection compo- 
nent with the other lower energy bands. In contrast to the previous 
findings, the CCFs of the reflection component (E4) with the lower 
continuum (El), and iron line region (E2) points at possible time 
shifts (< 1000s) as suggested in Fig. |l(]. Fig. |l4] b illustrates that the 
CCFs of a time shifted light curve is marked by a non-symmetric 
bodily shift of the CCF, whereas this is not the case of the CCF of 
a light curve with a phase shift. 

Fig. [n| shows that neither a time nor phase shift is seen in a 
comparison of the iron line region (E2) with the lower continuum 
(El), and the reflection component (E4) and upper continuum (E3). 
CCFs of the individual energy bands with itself for all energy bands 
mentioned thus far look nearly identical to Figs. 

We note that errors are such that we are unable to make defini- 
tive statements either about a phase or time lag at this time. 

5.1 Large scale bumps versus small scale flicker 

It is interesting to compare our time-lag results with those of Nowak 
& Chiang (1999) and Reynolds (1999). Nowak & Chiang use sim- 
ilar cross-correlation techniques to those employed here to search 
for time lags between the soft band (0.5-1 .0 keV) ASCA light curve 
and the hard band (8-15 keV) RXTE light curve. They found that 



the hard band lags the soft band by 1.6 ± 0.5 ksec. In this work, we 
find lags < 1000 s between RXTE bands El (2-4.5 keV) and E4 
( 10-20 keV). Noting that thermal Comptonization predicts a time 
lag that varies logarithmically with energy, our results are consis- 
tent with those of Nowak & Chiang. 

Reynolds (1999) uses an interpolation method to constrain 
trial transfer functions linking two given bands. He found a lag of 
50-100 s between the 2-4 keV and 8-15 keV RXTE bands, rather 
smaller than that found here. To reconcile these results, one must 
appreciate that these methods probe lags at different Fourier fre- 
quencies. The CCF methods tend to probe lags across a broad 
spectrum of Fourier frequencies. Due to the red nature of the 
power spectrum, such methods are naturally biased towards the 
lower Fourier frequencies. The method of Reynolds (1999), in- 
stead, probes the higher Fourier frequencies since he uses fairly 
spiky trial transfer functions. Hence, to paraphrase these technical 
results, the rapid flickering seems to get transmitted up the observed 
energy spectrum with a smaller time lag (by an order of magnitude) 
than experienced by the slower variations. 



6 POWER SPECTRA AND PERIODICITY ? 

While spectral studies of X-ray variability in time sequence may 
hold the key to understanding the underlying processes that are re- 
sponsible for producing the observed dramatic flux changes, it is 
insufficient for constraining the size of the emitting region in the 
absence of a good understanding of the flare mechanisms. The line 
profile obtained from ASCA observations suggest that the X-ray 
emission originates from ~ 10-20 gravitational radii of the black 
hole. This together with a periodic signal can constrain the size of 
the emitting region. Alternatively, we can attempt to estimate the 
black hole mass by assessing where the break frequency in the PDS 
occurs (assuming that the break frequency scales with mass). 

We calculate the power density spectrum using the Lomb- 
Scargle method (Lomb 1976; Scargle 1982; Press et al. 1992) ap- 
propriate for unevenly sampled data, and fit a power law slope in- 
dependently to the RXTE and ASCA data between 10~ 5 and 10~ 4 
Hz. (We ignore data above the latter in order to avoid contamina- 
tion to the fit from the 96 minute orbital period of RXTE .) Figs. [l8| 
show that respectively, /~ 1,3 and / _1 ' 5 is sufficiently representa- 
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Figure 15. Data points indicate possible phase lag in (a) CCF between the El (2—4.5 keV) and E3 (7.5—10 keV) light curves and (b) CCF between the E2 
(5—7 keV) and E3 (7.5—10 keV) light curves; superimposed on each CCF is its mirror image as represented by the solid lines. For both cases, the pattern 
between the CCF and its mirror image is similar to that of Fig. 14p for the simulated data, when a phase shift of <f> = 0.6 rad is applied. 
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Figure 16. Data points indicate possible time lag in (a) CCF between the El (2—4.5 keV) and E4 (10—20 keV) light curves and (b) CCF between the E2 
(5—7 keV) and E4 (10—20 keV) light curves; superimposed on each CCF is its mirror image as represented by the solid lines. For both cases, the pattern 
between the CCF and its mirror image is similar to that of Fig. |l4p for the simulated data, when a time shift of < 1000s is applied. 
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Figure 17. There are no indications for obvious shifts, either in phase or time for the CCF of the (a) El (2—4.5 keV) and E2 (5—7 keV) light curves, and (b) 
the E3 (7.5-10 keV) and E4 (10-20 keV) light curves. 



tive of the RXTE and ASCA data down to ~ 4 - 5 x 10" 6 Hz, 
where the break in the power spectrum may occur. (We note that 
this value is only given as a limit to fb r - while there appears to 
be no additional evidence for a break below 4 — 5 x 10~ 6 Hz, the 
observations are insufficiently long to claim a definitive determina- 
tion.) This is consistent with the findings of Hayashida et al. (1998) 
for MCG— 6-30-15. Since the count rate throughout these observa- 
tions remain steadily between 6-24 ct s _1 for RXTE and 0.5-3 ct 
s _1 for ASCA , the findings above would suggest that large scale 
power does not exist in abundance even though much shorter time 
scale variability is highly evident. 

We note that Papadakis & Lawrence (1993) caution against 
standard Fourier analysis techniques for estimating power spectra, 
in the form of a bias of the periodogram due to a windowing effect, 
in addition to a possible 'red noise leak'. The former is tied to a 
concern that the sampling window function can alias power from 
frequencies above and below the central frequency ui p , thereby dis- 
torting the true shape of the power spectrum; the latter 'red noise 
leak' effect is concerned with a transfer of power from low to high 
frequencies. (The problem of the 'red noise leak' was first noted 
by Deeter & Boyton (1982) and Deeter (1984).) However, we con- 
clude according to subsequent equations (6) and (7) that the effects 
of this bias is negligible for the data sets in question. The expected 
value of the periodogram is defined such that : 



RXTE and ASCA light curves 



E[I(u 



f(oj)Fpf(LL> — uj p )dw 



(5) 



where Fn(u — oj p ) is the Fejer kernel (Priestley 1981) which as- 
sumes the shape of the sine 2 function for frequencies ~ ujp.lt 
follows that the bias is : 

6(w) = E[I(up)] - f(u) (6) 

For lengthy time series (e.g. 400 ks), the mean value of the pe- 
riodogram would tend increasingly more towards the true value of 
the power spectrum at frequency oj p , as the Fejer kernel becomes 
increasingly concentrated around this frequency. In any case, we 
only wish to note that a potential break is seen in the power spec- 
trum at 4 — 5 x 10~ 6 Hz. (A similar break is noted in this object 
from Ginga data by Hayashida et al. 1998.) 

We wish additionally to point out an interesting possibility for 
a 33hr periodicity. This is illustrated in Fig. with 33hr interval 




Figure 19.RXTE (upper) and ASCA (lower) lightcurves of MCG-6-30-15 
from August 1997. Numbering the six tickmarks from left to right, ticks 2 to 
5 correspond to small sharp flares seen in the RXTE data; 2, 4 and 5 are also 
seen in the ASCA data. There is a peak near tick 6 and tick 1 corresponds 
to the (jagged) exit from a dip. Ticks 2 to 5 are regular to within about 1 per 
cent. 



tickmarks superimposed upon the ASCA and RXTE light curves. 
This has been determined by taking the mean of the times of the 2 
brightest flares in the RXTE light curve (the peak of the RXTE flare 
XI and i5, shown in Fig. 1). We note that the power of these peaks 
is not sufficient to be significant in Lomb-Scargle power spectra 
(i.e. the sharp peaks do not carry very much power). Accordingly, 
we have not attempted to quantify the significance of the peaks, and 
in part also because of the red noise nature of the PDS. We merely 
point out that 5 out of the 6 tickmarks in the RXTE light curve with 
flux > 14 ct s -1 occur in 33 hr intervals, with the same trend seen 
in the ASCA light curve. 



7 DISCUSSION 

There is evidence that the observed spectral variability is complex. 
However, all evidence points predominantly to a steepening of the 
spectral index with increasing flux. (Other AGNs that have exhib- 
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Figure 18. Lomb power as function of frequency for (a) RXTE 2-60 keV PCA energy band and (b) ASCA 0.6-10 keV. A power law model was used to fit the 
data; best fit values for the slope are labeled in the figures. 



ited spectral variability include NGC7314 Turner 1987; NGC2992 
Turner & Pounds 1989; NGC4051 Matsuoka et al. 1990; NGC3227 
Turner & Pounds 1989, Pounds et al. 1990; 3C273 Turner et al. 
1989, 1990; NGC5548 Nandra et al. 1991, Chiang et al. 1999; 
NGC4151 Perola et al. 1986, Yaqoob & Warwick 1991;) 1H0419- 
577 Guainazzi et al. 1998.) This may be indicative of changes in 
the temperature or optical depth of the Comptonizing medium or 
of the soft local radiation field. For instance, we can postulate that 
during periods of intense flux, a substantial amount of the hard X- 
rays are absorbed by the disk, and thermalized resulting in a source 
of soft photons. These will pass through the corona and Compton 
cool it, thereby giving rise to a steeper spectral slope. However, 
unless we fully understand the nature of the competition between 
coronal heating and Compton cooling, it is not clear what physics 
dominates to give the observed behaviour. 

The ejection model of Beloborodov (1998) can explain the 
observed relationship between reflection fraction and spectral index 
shown in Table [| In this model, observed spectral features are due 
to a non static corona, in which flares are accompanied by plasma 
ejection from the active regions. In other words, the bulk velocity 
( P = v/c ) in the flare is > 0. As the 'blobs' move away (are 
ejected) from the disk, a decrease in reflection is seen as a result of 
diminished reprocessing due to special relativistic beaming of the 
primary continuum radiation away from the disk. The model does 
not however account for the constancy of the iron line flux. 

We note that the possibility that some of the the iron line and 
reflection may be due to some distant material such as e.g. a torus is 
ruled out for this data based on (1) 196, and 199 findings for a non- 
constant, and sometimes absent narrow core seen in long ASCA 
observations of MCG-6-30-15 in 1994 and 1996, and (2) RXTE 
findings in this paper that the reflection fraction is seen to increase 
with flux. 



7.1 Reflection and the iron line 

A major result from the present observations is the enigmatic be- 
haviour of the iron line; the inverse proportionality between Wkc 
and R and of Rnorm with Fn a , GF91 predict that the bulk of the 
line arises from fluorescence in optically thick material. In this sim- 
ple reflection picture, we would expect R and Wko to be propor- 
tional to each other provided (1) the Compton reflection continuum 
does not dominate the iron line region, (2) the state of the illumi- 



nated region does not change, and (3) the primary continuum has a 
fixed spectral shape. We note for the last point that GF91 point out 
that differences in the photon index up to Ar ~ 0.2 (as compared 
to our flux-correlated findings for Ar ~ 0.06) will contribute less 
than a 10 per cent effect to Wkc- 

This lack of proportionality between Wkc and R (and 
Rnorm with Fko), in conjunction with an apparently constant iron 
line may point at changes to the ionization of the disk in MCG— 6- 
30-15, which is explored further below. We note that tentative evi- 
dence for observed changes in Fkc does exist during time inter- 
vals surrounding flare events. Such variability is clear from the 
ASCA analysis (199). It is possible that F^a changes on time 
scales shorter than is resolvable from time-averaged spectra; in 
other words, changes to Fn a may only become resolvable during 
bright flares. 

We note that Chiang et al. (1999) find similar results for the 
constancy of the iron line and inverse proportionality between Wk q 
and R in their multi-wavelength campaign of the Seyfert 1 galaxy 
NGC5548. 



7.2 A simple model for the observed spectral variations 

We propose the following model in order to explain some of the 
enigmatic properties of the observed variability phenomena. Spec- 
tral variability is no doubt complex, and does not conform to the 
present picture of a cold disk geometry for MCG— 6-30-15. If how- 
ever the variable emission from MCG— 6-30-15 is from a part of 
the disc which is more ionized, say with an ionization parameter 
£ ~ 100 (see spectra in Ross & Fabian 1993, Ross et al 1999), then 
the reflection continuum will respond to the flux while the iron line 
does so only weakly. At that ionization parameter the iron line can 
be resonantly scattered by the matter in the surface of the disc and 
its energy lost to the Auger process (Ross, Fabian & Brandt 1996). 

The more highly ionized region could either be the innermost 
regions of the disc, perhaps within say 6r g , or the regions directly 
beneath the most energetic flares. We note that flux-correlated 
changes in the surface density of the disk can lead to changes to 
ionization states without incorporating large changes to luminosi- 
ties (Young et al. 1999, in preparation). 

To illustrate in the context of the RXTE light curve shown in 
Fig. [TJ, assume that the flux below 10 cts -1 reflect physical pro- 
cesses that occur within the radius 6 
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is enhanced in the > lOcts -1 variability, by flares within 6r g , 
where the Auger destruction effect becomes important. Accord- 
ingly, this will lead to observable changes in reflection (i.e stronger 
reflection during the higher flux periods), with minimal changes in 

FKa- 

We note that our interpretation that variability largely comes 
from within the innermost stable orbit for a Schwarzschild black 
hole may be consistent with the scenario for a very active corona 
(and hence strong hard X-ray emission) within 6r 9 , proposed by 
Krolik (1999). In this model, magnetic fields within the radius of 
marginal stability are strong and amplified through shearing of their 
footpoints, which can enhance variability. 

7.3 Implications for Mass Estimates 

7. 3. 1 Constraints from spectral studies 

The constancy of the iron line on day-to-day scales suggests that 
the timescale for variability (i.e. the observed periods for which 
dramatic flux changes are observed) we are naively probing are 
much larger (in the 'standard' scenario) than the fluorescing re- 
gion. In other words, slow changes would imply on a naive model 
much larger crossing times and hence large regions for the crossing 
times of the continuum. A light-crossing time of the fluorescing re- 
gion larger than ~ 50 ks (assuming an average radius ~ 20r s ) will 
lead to an estimate for the black hole mass ~ 1O 8 M . Reynolds 
(1999) points out however that the bulge/hole mass relationship of 
Magorrian et al. (1998) implies a much lower mass estimate for 
MCG-6-30-15, by an order of magnitude, of about ~ 10 7 M Q . 

In the scenario of the simple model presented above, and given 
evidence for short timescale variability of the iron line (here and 
199) as well as the location of the flare line found in 199, the con- 
stancy of the line suggests that the timescale for variability for 
which we are probing is much smaller than the fluorescing region 
and reconciles the above mass problem. 

7.3.2 PDS : Analogies with Galactic Black Hole Candidates 

Of further interest is the apparent break in the power spectrum of 
MCG-6-30-15 seen in both the RXTE and ASCA data. The ori- 
gin of the break is not yet known (but see e.g. Edelson & Nandra 
1999; Poutanen & Fabian 1999; Kazanas, Hua, & Titarchuk 1997; 
and Cui et al. 1997 for possible explanations) , but does provide 
a useful means to determine the black hole mass, through scaling 
from similar breaks in the power spectrum of the famous galactic 
black hole candidate (GBH) Cygnus X-l, and other objects like it. 

The behaviour of the PDS in MCG-6-30-15 is not unlike that 
of GBHs in the 'low' (hard) state (see e.g. Belloni & Hasinger 
1990; Miyamoto et al. 1992; van der Klis 1995). Power law slopes 
(with form /") of order a ~ -1 to -2 are observed at high frequen- 
cies and flatten to ~ at lower frequencies. If we bridge the gap be- 
tween AGNs and GBHs and assume that similar physics are at play, 
we can make predictions for the black hole mass in MCG— 6-30-15 
(using the values for the cutoff frequencies ft r ) by a simple scaling 
relation with Cygnus X-l. Belloni & Hasinger (1990) report that 
fbr ~ 0.04-0.4 Hz for Cygnus X-l; for MCG-6-30-15, we find 
evidence that f br ~ 4 — 5 x 10~ 6 Hz. The resulting ratio between 
the 2 cutoff frequencies is ~ 10 4 — 10 5 . Herrero et al. (1995) argue 
that the black hole mass in Cygnus X-l is ~ 1OM , which leads 
us to conclude that the mass of the black hole in MCG— 6-30-15 is 
~ 10 5 — 10 6 , smaller than anticipated. 

Our mass estimate agrees with that of Hayashida et al. (1998) 



and Nowak & Chiang (1999) who also used the break frequency 
and scaling arguments. However, such mass estimates should be 
treated with extreme caution. The break frequencies in any one 
given GBHC can vary by one or two orders of magnitude depend- 
ing upon the exact flux/spectral state of the source. Given that we 
do not know how to map AGN spectral states into analogous GBHC 
states, the mass estimate derived above (and that of Nowak & Chi- 
ang 1999) will also be uncertain by up to two orders of magnitude. 
Additionally, it is not entirely clear what timescales to identify the 
cutoff frequency with. A small black hole mass (i.e. < 2 x 1O 6 M ) 
would also imply the presence of a super-Eddington black hole in 
MCG-6-30-15. 

7.3.3 A possible 33 hour period 

Finally, we address what implications a 33 hr period would have on 
the black hole mass in MCG— 6-30-15. If we make the assumption 
that this is the orbital time scale for e.g. a flare to circumnavigate 
the black hole in MCG— 6-30-15, then we can estimate the mass 
via the relation : 

M BH (R) = 3.1 x 10 7 t orb (-A-)- 15 M (7) 

LUH S 

where R is the distance from the center, and R a = 2r 9 is the 
Schwarzschild radius (the gravitational radius of the black hole 
r g = GM/c 2 , and t or b is in days). 

The diskline model constrains Ri n and R ou t assuming some 
power law emissivity function (oc R~ a ) that declines to larger 
radii. (We note that beyond r out the line emission is negligible.) 
Accordingly, we expect that most of the power is concentrated in 
the inner radii. 196 and 199 constrain using time averaged ASCA 
data r in ~ (6.7 ± 0.8)r 9 for MCG-6-30-15. This combined 
with t or b ~ 33hr (= 1.375 days) give a mass for the black hole 
in MCG-6-30-15 ~ 2.6 x 10 7 M Q . 



8 CONCLUSION 

We summarize below the spectral and timing results of this pa- 
per. It is clear that complicated processes are present, the nature of 
which is not obviously apparent, and may prove to be a challenge 
to present theoretical models. 

• Hardness ratios reveal that spectral variability may be largely 
attributed to changes in the intrinsic photon index. In particu- 
lar, spectral hardening is observed during periods of diminished 
intensity in comparisons of the (7.5—10 keV) upper continuum, 
(5—7 keV) iron line region, and (10—20 keV) reflection hump with 
the (3—4.5 keV) lower continuum. Particularly hard spectra are 
noted in a time interval corresponding to ~ 260 ks that begins 
shortly after the hard RXTE flare. 

• We find from flux correlated studies that changes to the photon 
index are evident. In particular, r3_io steepens while Tio-20 flat- 
tens with flux; for a doubling of the flux, Ar3_io ~ 0.06 and 
ATio-20 ~ 0.3. This coupled with findings for a constant iron line 
can contribute to the reduced fractional variability in the iron line 
band noted by Reynolds (1999). We note that changes to Tio-20 
are significant only with large changes in flux whereas changes in 
r3_io are apparent even with subtle changes in flux. This point is 
well illustrated in detailed studies of the time intervals surrounding 
the ASCA and RXTE flares. Nevertheless, it would appear that both 
changes in the intrinsic power law slope (reflected by changes to 
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r3-io), and reflection (reflected by changes to r 10-20) both con- 
tribute in varying degrees to the overall spectral variability. 

• We find curiously that the iron line flux is consistent with being 
constant over large time intervals on the order of days (but this is 
not the case on much shorter time intervals of order < 12 ks), and 
the equivalent width anitcorrelates with the continuum flux. (Ob- 
served changes to Fk q on short time intervals are summarized in 
the next point.) This may point at evidence for a partially ionized 
disk. 

• From concentrated studies of the time intervals surrounding the 
RXTE and ASCA flares, we find tentatively that Fko shows a no- 
ticeable increase after the flare events. (This is less significant for 
the periods of the RXTE flare.) This may be an indication that we 
are witnessing some type of response to the flare. We note that Fkc 
is comparably high (~ factor of 1 .7 larger) during the times sur- 
rounding the RXTE flare, and times following the ASCA flare, in 
contrast to time averaged analysis of flux-correlated data. 

• We find tentative evidence from cross correlation techniques 
for a possible phase lag comparable to <j> ~ 0.6 between the 
(7.5 — 10 keV) upper continuum, and (5—7 keV) iron line band and 
(3—4.5 keV) lower continuum. 

• CCFs further reveal possible time lags (time delays < 1 ks) be- 
tween the (10—20 keV) reflection hump and iron line band, and 
reflection hump with lower continuum. 

• We report an apparent break (from f° to ~ / _1 ' 5 ) of MCG— 6- 
30-15 at ~ 4 - 5 x 10~ 6 Hz (~ 56 - 70 hrs) seen by both ASCA 
and RXTE . Scaling with the mass of the GBH Cygnus X-l gives 
a smaller than expected black hole mass of - 10 5 - 1O 6 M for 
MCG— 6-30-15. However, this is unlikely to be a proper estimate 
of the mass for the black hole in MCG— 6-30-15. (A black hole 
mass < 2 x 1O 6 M would make the black hole in MCG-6-30-15 
super-Eddington.) 

• We report on the possibility for a 33 hr period seen in both the 
ASCA and RXTE light curves. This combined with a value of ~ 
7r g for the inner radius, implies black hole mass ~ 2.6 x 10 7 M Q 
for MCG-6-30-15. 
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